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> 

in 

in ■ We reexamine contributions of higher if-resonances to the radiative rare 

00 ' decays b — > 57 in the limit where both 6- and s-quark are considered heavy. Us- 



ing the non-relativistic quark model, and the form factor definitions consistent 
with the HQET covariant trace formalism, we find significant disagreement 
with previous work which also used heavy quark symmetry, and excellent 

■ agreement with experimental results. In particular, the two largest fractions 
of the inclusive b — > 57 branching ratio are found to be (16.8 ± 6.4)% for 
B A"* (892)7 and (6.2 ± 2.9)% for B K% (1430)7 decays. We also com- 
pare the contribution from the radiative decays into the eight iT-meson states 

■ to the inclusive experimental b — > S7 mass distribution. 



1 Introduction 



Flavor changing neutral current transitions involving the 5-meson provide a unique 
opportunity to study the electroweak theory in higher orders. Although transitions 
like b — > s 7 , b — > se + e~ , and b sg vanish at the tree level, they can be described 
by one loop ("penguin") diagrams, in which a W~ is emitted and reabsorbed 
These processes occur at a rate small enough to be sensitive to physics beyond 
the Standard Model f2|. Similar flavor violating processes in the iT-meson system 
have the disadvantage that non-perturbative long distance effects are quite large, 
and it is difficult to extract the quark level physics from well known processes like 
K + — > 7r + e + e~. 

Among all rare B decays, radiative processes B — > X s7 (especially decay B — > 
K*(892) r y) have received an increasing attention, because of the experimental mea- 
surement of the B — > K* (892)7 exclusive branching ratio ||, 

BR(B -> (892) 7 ) = (4.5 ± 1.5 ± 0.9) x 1(T 5 . (1) 

which has been recently updated to 

BR(B -> iT(892) T ) = (4.3l};J ± 0.6) x 10" 5 , (2) 
and also of the inclusive rate f| , 

BR(B -> A s7 ) = (2.32 ± 0.57 ± 0.35) x 10" 4 . (3) 

Several methods have been employed to predict exclusive B —>■ K*(892) / y decay rate: 
HQET QCD sum rules P1-|T3]1, quark models |]l4]|-pl|, bound state resonances 

p5| , and Lattice QCD PB|-[29|. The theoretical uncertainty, which was originally 
of two orders of magnitude, has been greatly reduced in the more recent studies. 
However, there is still a large spread between different results. 

In this paper we follow the approach of @, [Tj, in which both b- and s-quark are 
considered heavy. In the heavy quark limit the long distance effects are contained 
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within unknown form factors, whose precise definition consistent with the covariant 
trace formalism [3C]-[32] has been clarified only recently This is precisely the 



reason why our results substantially differ from 0, even though we use the same non- 
relativistic quark model for the wave functions of the light degrees of freedom (LDF). 
Our results show that the ratio of the exclusive B — > K**^ to the inclusive decay 
rate B — > X s 7 was underestimated for the channel B —>■ K*(892) / y ((16.8 ± 6.4)% 
as opposed to (3.5 — 12.2)% from J7|), and significantly overestimated for the decay 
B -> #5(1430)7 ((6.2±2.9)% as opposed to (17.3-37.1)% from 0). We emphasize 
that our prediction for the decay B — > K* (892)7 is in agreement with experimental 
result of (19±5)%. Although other exclusive decays have not yet been identified, we 
have compared with experiment the contribution from the eight B — > K**^ decays 
to the inclusive B — > X s 7 mass distribution. 

The paper is organized as follows: in Section we restate the theoretical frame- 
work for the B — > K** / y decays. Section contains a discussion of the form factor 
calculation. The expressions for the form factors given in are evaluated in terms 
of the wave functions and energies of the light degrees of freedom in the meson rest 
frame. We discuss here the model used in establishing the LDF wave functions and 
energies. An extensive literature exists in this subject, so we have attempted to set 
our results in context with previous calculations in Section [|. Our conclusions are 
summarized in Section 

2 Theory of B — > K**j decays 

The effective Hamiltonian for the decays B — > A s 7 can be found in many places, e.g. 



34j-|36fl. It is derived by integrating out the top quark and VT-boson at the same 
scale fi ~ Mw An appropriate operator basis for the effective Hamiltonian consists 
of four-quark operators and the magnetic moment type operators of dimension six 
(Oi — 0$). Higher dimensional operators are suppressed by powers of the masses of 
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the heavy particles. For the B — > K**^ decays only the operator 7 contributes, so 
that 

4G F . 



Heff = --j=-V tb V t :C 7 (m b )0 7 (m b ) . (4) 



Here, 7 is given by 



7 = ^F, v [m b sa^{l + l5 )b + m s sa^{\ - 75)6] , (5) 



with a^ u = I [7^,7^]. The explicit expression for the Wilson coefficient C 7 (m b ) 



-"li- 



as a function of -jjt- can be found in 36, 37 . The value of C 7 can be calculated 



perturbatively at the mass scale /1 = Mw- The evolution from M\y down to a mass 
scale ji = m b introduces large QCD corrections. This proceedure also introduces 
large theoretical uncertainties, primarily due to the choice of the renormalization 
scale /1 (taken above as m b ), which can be as large as 25% [ j36|| . 

As proposed in || 0, we evaluate the hadronic matrix element of 7 between 
a S-meson in the initial state, and a generic i^**-meson in the final state, in the 
heavy quark limit for the b- and s-quarks. Matrix elements of bilinear currents of two 
heavy quarks ( J(q) = Q'TQ) are most conveniently evaluated within the framework 
of the trace formalism, which was formulated in []30l |3l| and generalized to excited 



states in [32]. Denoting to = v ■ v', where v and v' are the four- velocities of the two 



mesons mesons, we have 



(y'{v')\J{q)\^{v)) = Tr[M'{v')TM{v)]M{uj) , (6) 

where M' and M denote matrices describing states and ^(v), M = 7°M"'~7 , 

and M.(oj) represents the LDF. For all transitions considered in this paper matrices 



M and M', as well as definitions for Ai(uj), can be found in 0, [33|]. Using (^j), we 
can write 

(K**j\O r (m b )\B) = ^ V ^Tr[M>(v')W»M(v)]M(uj) , (7) 
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where the factor q v = rriBV v — mK**v' v came from the derivative in the field strength 
, rj^ is the photon polarization vector, and 



of 



Q!»> = m B a^{\ + 75) + m K **a^{l - 75) • 



(8) 



Expression (^) can be further simplified using j6M(v) = M(v). 

Now, using the mass shell condition of the photon (q 2 = 0), and polarization 
sums for spin-1 and spin-2 particles, we obtain the following decay rates M: 



T(B - 


-> K* (892) 7 ) 


= 




2 i[(i-y) 3 (i + y) 5 (i + y 2 )], 
y 


(9) 


F(B - 


♦ ifi(1270)7) 


= ni^i 




\ 2l -[(l-y) 5 (l + y) 3 (l + y 2 )], 


(10) 


T(B - 


♦ tfi(1400)7) 


= m F [ 


U ) 




(11) 


T(B - 


>#*(1430) 7 ) 


= n\£ F ( 


u>) 


2 7^[(l-2/) 5 (l+2/) 7 (l+y 2 )] , 

8r 


(12) 


T(B - 


> #*(1680)7) 


= m G ( 


u) 


l 2 ^[(i-y) 7 (i + l/) 5 (i+2/ 2 )]> 


(13) 


T(B - 


* # 2 (1580)7) 




u) 


\ 2 ^- 3 [(i-y) 7 (i + yf(i + y 2 )} , 


(14) 


T(B - 


> A"* (1410)7) 


= m C2 


(u 


)\ 2l -[(l-y) 3 (l + y) 5 (l + y 2 )}, 


(15) 


T(B - 


* #i(1650)7) 


= n\u 


(CO 


)| 2 I[(l_ 2/ )5 (l + y) 3 (l + y 2 )] j 


(16) 



where we used abbreviations 

y 
n 



TXlK** 

m B 
a 

128tt z 



G 2 F ml\V tb \' 2 \V ts \ 2 \C 7 (m b )\' 



(17) 
(18) 



and the argument of the Isgur-Wise (IW) functions is fixed by the mass shell con 
dition of the photon (q 2 = 0), 

1 + y 2 



uJ 



2y 



(19) 



Note that in the expressions for the decay rates (P)-([TB[) given in 0, a factor of 
(1 — y 2 ) was omitted. Also, as observed in 0, since decays into the states belonging 
to the same spin symmetry doublet are described by the same Isgur-Wise function, 
and since in the heavy-quark limit the two members of a spin doublet are degenerate 
in mass, from (§)-(|16|) one has 

T(B #*(1430)7) « 3r(S if i(1400)7) , (20) 
r(5 ^ iT 2 (1580)7) ~ 3r(5 -> ifi (1680)7) • (21) 

As indicated, these relations are only approximate due to a large breaking of the 
spin symmetry for the s-quark. 



3 Model for the Isgur-Wise functions 

As already mentioned, even though we use the same non-relativistic quark model, 
our calculation differs significantly from |7j in evaluation of the IW form factors 
needed for the decay rates. Assuming that we can describe heavy-light mesons 
using a non-relativistic potential model, the rest frame LDF wave functions (with 
angular momentum j and its projection Xj), can be written as 

4>j\? ( x ) = X! R aL{r)Y LmL (Q)x ms ( L > m L, ^,m s \j, XfL, -) , (22) 



m L ,m s 



where Xm s represent the rest frame spinors normalized to one, xL>Xm a = &m' s ,m 3 , 
and a represents all other quantum numbers. According to instead of the 



simple overlap of the two wave functions, the form factor definitions should include 
a Lorentz invariant factor in front of the overlap of the two wave functions describing 
the initial and the final states of the LDF. Also, following the suggestion of [|38[| , 
overlaps of the two LDF wave functions can be done in the Breit frame (v = — v'), 
where the boost factors (connecting the moving to rest LDF states) cancel out. All 
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this leads to the following expressions [33] valid for the non-relativistic quark modelQ 
(suppressing quantum numbers a' and a, and using the notation of 0): 

CcM = -?—(j (ar))ao , 07 -(07, IT), (23) 

LO + 1 2 2 2 

£ E (u) = J_ - 0i(or))i , 0T-(0+ 1+), (24) 

VW — 1 2 2 2 

M^) = J / . 07-(lt,2t), (25) 

V or — 1 uj + 1 2 22 



i G {u) = ^-(j 2 {ar)) w , 07 -(17,27), (26) 



where (denoting the energy of the LDF as E t 



a = {E 9 + E f 9 ) ) j^ Tl , (27) 

and 

(Hr))il = J r 2 drR* a , L ,(r)R aL (r)F(r) . (28) 

Note that (p3|)-(f26|) include transitions from the ground state into radially excited 
states. If the two 3 = \ states are the same, E'= = Eg and £c7 is normalized to one. 
The above expressions should be compared with the ones used in [0 (putting a tilde 
over the form factors to avoid confusion), 

ic(cu) = O'o(ar)) o, (29) 

i E (u) = >/30i(ar)) 10 , (30) 

£ F ( U ) = y/l(j x {ar)) 10 , (31) 

ia(u) = V5(j 2 (ar)) 20 , (32) 



with the definition 



E'VJ^I . (33) 



1 As pointed out in p3| , models based on the Dirac equation with a central potential lead to 
the same expressions for the IW functions. 
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For the numerical estimates we employ the model used in [pj (usually referred 
to as the ISGW model), the Schrodinger equation with 

V(r) = -— - + c + br . (34) 
3r 

With sensible choice of parameters, this simple model gives quite reasonable spin- 
averaged spectra of bd and sd mesons up to L = 2. However, instead of just using 
a single harmonic oscillator wave function (as was done in 0), for the radial wave 
function of the LDF, we numerically solve the Schrodinger equation. To determine 
the parameters of the model, we fix b = 0.18 GeV 2 (which was also used in P9j), 
and vary a s and c for a given value of m u ^ (in the range 0.30 — 0.35 GeV), and 
m s (in the range 0.5 — 0.6 GeV), until a good description of the spin averaged 
spectra of i^-meson states is obtained. Following this proceedure, our a s ranges 
from 0.37 to 0.48, while c takes values from —0.83 GeV to —0.90 GeV. These 



parameters are in good agreement with the original ISGW values (a s = 0.50 
and c = —0.84 GeV for m U) d = 0.33 GeV and m s = 0.55). We emphasize that the 
original ISGW parameters give results that are well inside the ranges for all decays 
quoted in this paper. By varying the c- and 6-quark masses we could also obtain 
good spin averaged description of the B and D mesons. However, to be consistent 
with heavy quark symmetry, the wave function for the B meson was chosen to be 
the same as the one obtained for the spin averaged (ground state for L = 0) K and 
if* (892) mesons. 

To completely define our proceedure, we have to specify how the LDF energy Eq 
was determined. In []7| for a given i^**-meson the LDF energy was defined as 

= rnK~*m^ 
m s + m Utd 

This definition was proposed to account for the fact that s mesons aren't particularly 
heavy. On the other hand, a definition that is consistent with heavy quark symmetry 
is 

Eg = m K ** - m s . (36) 
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It should be noted that these two expressions are not equivalent in the heavy quark 
limit. In order to explore the sensitivity of our results on the choice of Eq, we have 
repeated all calculations employing both of these two definitions, and in the final 
results we have quoted the broadest possible range obtained for the form factors 
(and for all other results). Finally, Eq for the B meson has been taken to be the 
same as Eq for the i^*(892) meson, consistent with heavy quark symmetry. It turns 
out that this is actually a very reasonable assumption. The range of Eq that was 
used here for B and K*(892) meson was from 0.296 GeV to 0.396 GeV. On the 
other hand, from the CLEO data on the semileptonic B decays p0| , and the LQCD 



heavy-light wave function [4T|, it was estimated |^2[ that in B systems Eq ranges 



from 0.266 GeV to 0.346 GeV. 

We believe that the proceedure outlined above enables us to estimate a reasonable 
range for the unknown IW form factors in a physically more acceptable way than 
it was done in pj, by simply varying the scale parameter of the single harmonic 
oscillator wave function. 

4 Our results and comparison with previous in- 
vestigations 

In Table [I] we present our results for the range of (absolute) values of the form factors 
at the indicated value of lu, for the ratio R = p^z^ ^ ? an d for the branching 
ratio BR(B — > K**j), for the various i^**-mesons. The inclusive branching ratio 
B — > X s 7 is usually taken to be QCD improved quark decay rate for b — *> 57, which 
can be written as [El, El 



1 1 

Ttl Ttl 

T(B -> X sl ) = 40(1 - ^) 3 (1 + ^4) . (37) 



The leading log prediction for BR(b — > s-y) is (2.8 ± 0.8) x 10 4 (pj| B_7|], where the 
uncertainty is due to the choice of the QCD scale. The next-to-leading order terms 
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that have been calculated tend to reduce the prediction to about 1.9 x 10 4 [p7[ . 
Both of these predictions are in excellent agreement with the recent experimental 
result of BR(b -> sj) = (2.32 ± 0.57 ± 0.35) x 10~ 4 §. For the numerical values of 
the B — > K**'-/ branching ratios given in Table p] we used the leading log result of 
BR(b -> S7) = 2.8 x 10" 4 . 

In order to make comparison of our results with previous calculations easier, we 
have tabulated our results together with results of |7]] and |IB| in Table As far 
as we know, these two papers are the only ones that have dealt with radiative rare 
B decays into higher X-resonances. There has been much more work done on the 
decay B — > K* (892)7), and we have tabulated some of these results in Table |3|. As 
one can see from Table ^ the predictions for this particular ratio ranges from a 0.7% 
fl22| to 97.0% 0. The data suggest a value of (19 ± 5)%. Note that our result of 
(16.8 ± 6.4)% is consistent with the data, unlike the values quoted in and [ |i~6[ . 
As far as decays into higher K resonances are concerned, our results are in general 
in much better agreement with |T3 than with M. In particular, the authors of [|7| 
emphasized a large branching ratio for the decay B — > .K^ (1430)7 ((17.3 — 37.1)%), 
while our results indicate a 3-6 times smaller value of (6.2 ± 2.9)%, a result which 
agrees with the one quoted in Jl6| (6.0%). Also note that our numerical results from 
Table ^ support relations (p0|) and (j2TJ). 

With the exception of the K* (892)7 channel, no other exclusive radiative pro- 
cesses have been identified so far. The inclusive radiative B — > Xg'j mass distribution 
has however been measured by CLEO ||, and is shown in Fig. |I[ We have normal- 
ized experimental data so that the integrated distribution gives unity. The i^*(892) 
peak is evident, but the higher mass contribution are not resolved. We have at- 
tempted to model this inclusive distribution by considering the contributions from 
each of the exclusive K**^ channels considered in this paper (and given in Table |2|). 

In order to compare our result to experiment, we replace a given Rj^** by a mass 
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distribution reflecting the finite total width Tk** of the K** resonance |[43||, 



dR{m Xs ) = Rk- T K ,*/2 

dm Xs £j tt (m Xs -m K ») 2 + (r^"/2) 2 ' 1 J 

The integrated distribution gives 

/^ m - = g B - (39) 

In Figure |l| we show the total resonance contribution (solid line) compared to the 
experimental inclusive B —>■ X s 7 mass distribution. The area of the resonance curve 
is 37.4% of the total inclusive rate (see Tables [I] or H). We see the general shape 
is correct, but it is difficult to make more quantitative statements due to the large 
errors involved. 



5 Conclusion 

In this paper we have reexamined predictions of heavy quark symmetry for the 
radiative rare decays of B-mesons into higher A-resonances. An earlier calculation 
@ suggested a substantial fraction ((17.3—37.1)%) of the inclusive b — > branching 
ratio going into the A;? (1430) channel, and only (3.5 — 12.2)% going into the A*(892) 
channel. Even though we used the same non-relativistic quark model, our calculation 
yields fractions of (16.8±6.4)% and (6.2±2.9)% for A^*(892) and A 2 *(1430) channels, 
respectively. Note that experimental results favor the value of (19 ± 5)% for the 
A* (892) channel. Besides a more careful treatment of the uncertainty in the wave 
functions of the light degrees of freedom, our calculation differs from in employing 
form factor definitions that are consistent with the HQET covariant trace formalism 
331. As a consequence of that, our results for all decay channels significantly differ 



from 0. The contribution of the eight A**7 channels to the inclusive B — > X s ~f 
mass distribution was compared with experiment. We find the general shape of the 



11 



mass spectrum to be correct, but due to the large errors involved one cannot reach 
more quantitative conclusions. 
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TABLES 



Table 1: Our results for the range of absolute values of the form factors at indicated 
value of u, for the ratio R = wgz^ ; an d f° r the branching ratio BR(B — > K**^), 
for the various i^**-mesons. For the calculation of branching ratios we used the value 
-> X s7 ) = 2.8 x 10" 4 [BH. 
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total 
















37.4 ± 13.6 


10.44 ± 3.78 
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Table 2: Comparison of our results for the ratio R = r^r^ $ with the previous 
work done in |7]] and [pT6 1 . Note that in the quark model calculations decay into 
the 1 Pi state is forbidden, because Oj is a spin-flip operator, and _fTi(1270) and 
-fTi(1400) are mixtures of 1 P\ and 3 Pi states. In [16| 3 Pi state had R = 6%. 



Meson 


State 


J p 


P[%] (this work) 


R[%] (ref. |7 


1) R[%] (ref. 0) 


K 


C 


0- 




forbidden 




if* (892) 


c* 


\~ 


16.8 ±6.4 


3.5 - 12.2 


4.5 


K*(1430) 


E 


0+ 




forbidden 




#i(1270) 


E* 


1+ 


4.3 ± 1.6 


4.5 - 10.1 


forbidden/6.0 


Kx(1400) 


F 


1+ 


2.1 ±0.9 


6.0 - 13.0 


forbidden/6.0 


K 2 *(1430) 


p* 


2+ 


6.2 ±2.9 


17.3-37.1 


6.0 


7^(1680) 


G 


V 


0.5 ±0.2 


1.0 - 1.5 


0.9 


K 2 (1580) 


G* 


2" 


1.7 ±0.4 


4.5 - 6.4 


4.4 


X(1460) 


c 2 


0- 




forbidden 




7^(1410) 


c; 


1- 


4.1 ±0.6 


7.2 - 10.6 


7.3 


K*(1950) 


E2 


0+ 




forbidden 




fTi(1650) 


E* 


1+ 


1.7 ±0.6 


not given 


not given 


total 




37.4 ± 13.6 


44.1 - 90.9 


29.1 
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Table 3: Comparison of our results for the ratio R = r | g ^ x ^ with several previous 
calculations. 



Author (s) 


Ref. 


R[%) 


Donnell (1986) 


[I4| 


97.0 


Deshpande et al. (1988) 


[I5| 


6.0 


Dominguez et al. (1988) 




I 


28.0 ± 11.0 


Altomari (1988) 


[I6| 


4.5 


Deshpande et al. (1989) 


[I9| 


6.0-14.0 


Aliev et al. (1990) 




i 


39.0 


Ali et al. (1991) 




i 


28.0-40.0 


Du et al. (1992) 






69.0 


Faustov et al. (1992) 


m 


6.5 


El-Hassan et al. (1992) 






0.7-12.0 


O'Donnell et al. (1993) 


[I7| 


10.0 


Colangelo et al. (1993) 


[I0| 


17.0 ± 5.0 


Ali et al. (1993) 






3.5 ± 12.2 


Ali et al. (1993) 




12 


16.0 ± 5.0 


Ali et al. (1993) 


[I8| 


13.0 ±3.0 


Ball (1994) 


[HI 


20.0 ± 6.0 


Nanson (1994) 


|13| 


16.0 ± 4.0 


Holdom et al. (1994) 




m 


17.0 ±4.0 


Atwood et al. (1994) 




m 


1.6-2.5 


Bernard et al. (1994) 




m 


6.0 ± 1.2 ±3.4 


Ciuchini et al. (1994) 




m 


23.0 ±9.0 


Bowler et al. (1994) 




m 


9.0 ± 3.0 ± 1.0 


Burford et al. (1995) 




m 


15.0-35.0 


Tang et al. (1995) 




N 


10.0-12.0 


this work 






16.8 ±6.4 
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FIGURES 



Figure 1: The experimental inclusive B — > X s ^ mass distribution measured at 
CLEO 0. The data have been normalized to unity. The curve is the sum of the 
exclusive K**^ channels from Table [I] as calculated by (|38|). 
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Figure 1 
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